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Abstract 

We compute the perturbative renormalization factors required to match to the continuum Isgur-Wise 
function, calculated using lattice Heavy Quark Effective Theory. The velocity, mass, wavefunction and 
current renormalizations are calculated for both the forward difference and backward difference actions 
for a variety of velocities. Subtleties are clarified regarding tadpole improvement, regulating divergences, 
and variations of techniques used in these renormalizations. 
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I Introduction 



The unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix is regarded as a crucial test of the standard 
model |, §; the precise determination of these matrix elements has received extensive experimental and 
theoretical scrutiny. The V c b CKM matrix element can be extracted from the reaction B — > D* IVi, if the 
theoretical factors in the decay rate can be reliably computed. The heavy quark effective field theory (HQET) 
formalism is well-suited to the analysis of this decay. The differential decay rate of the above process is 



-T(B^D*lV l ) = —^k(m B ,m D ,vv')\V cb \ z e(vv') (1) 

Cl(V ■ V j 407T 

where £(v ■ v') is a universal form factor, the Isgur-Wise function. The function k can be calculated in 
perturbation theory using various approximations [Q, [|. The Isgur-Wise function is a QCD matrix element 
that must be computed non-perturbatively. Previously and in a companion paper ^ pi, we discuss the 
numerical calculation of the Isgur-Wise function using lattice HQET. In this paper, we discuss the pertur- 
bative matching of lattice HQET to continuum HQET, which allows the conversion of the results from the 
numerical simulations into physical predictions. Specifically, we shall be matching from the lattice to the 
continuum matrix element, 



(D, v \J b ^ c (Q)\ B, v 1 ) = y/M D M B (»„ + t£) £(« • v') (2) 
where v and v' denote the 4-velocities of the c and b quarks, and 

J b ^ c {x)=c{x) lix (l-j 5 )b(x) (3) 

is the weak current for the transition of a bottom to a charm quark Q . 

The Isgur-Wise form factor describes the response of the quark-gluon sea surrounding the heavy quark 
due to a sudden change in velocity of the heavy quark when it decays. In HQET, the Isgur-Wise function is 
non-perturbatively equal to one at the point of zero recoil, v = v'; HQET does not constrain the Isgur-Wise 
function at non-zero recoil. Continuum perturbative corrections are required to obtain the zero recoil result 
in QCD; however, these are known to 2 loops 0. Unfortunately, there is no experimental data at zero recoil, 
so the experimental data is extrapolated Q to zero recoil in order to estimate V c b using Eq. (Q) . Knowledge 
of the functional form of the Isgur-Wise function would greatly aid this extrapolation. The Isgur-Wise 
function can be calculated non-perturbatively, in principle, from QCD for arbitrary recoil. In our companion 
paper [El, we describe our simulations that use lattice HQET to calculate the Isgur-Wise function. 

There have been previous calculations of the renormalization factors for lattice HQET. Unfortunately, 
not all of the perturbative factors required for our numerical simulations were calculated. After lattice HQET 
was introduced by Mandula and Ogilvie || , there were a number of concerns about the validity of the lattice 
HQET formalism Jj| |l(j, [llj. The consistency of lattice HQET was finally demonstrated by Aglietti 
in perturbation theory. However, Aglietti used a form of lattice HQET action that is less convenient for 
numerical simulation than the one originally used by Mandula and Ogilvie. The difference between the 



HQET actions was in the use of a forward or backward finite difference in the time direction (see Sec. til). 



Also, Aglietti considered only a special kinematic limit with one quark at rest and the other quark at finite 
velocity. Mandula and Ogilvie Jl3|] limited their work to the velocity renormalization factors for the forward- 
difference action (which we used in our simulations); they calculated neither the vertex function nor the 
wavefunction renormalization which are required to renormalize the lattice data. 

In this paper, we calculate the perturbative factors required to renormalize the Isgur-Wise function 
obtained from a lattice HQET simulation. The calculation includes two HQET actions: one with the 
forward time derivative and one with the backward time derivative. We follow the formalism developed 
by Aglietti Jl2| , but generalize Aglietti's expression for the vertex function to arbitrary input and output 
velocities (as is required for the analysis of the simulation data). We also include the effects of tadpole 
improvement and discuss a subtlety in the calculation of the vertex function. 
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Section [il] will provide a sparse review of continuum HQET in order to put the lattice calculation into 
context. Section III will describe the details of the velocity, mass, wavefunction, and vertex renormalizations 



for the lattice actions, including a discussion of the "reduced" results and an evaluation at non-zero recoil. 
Section IV will describe how these are combined into a single renormalization for the lattice current to be 
matched to the continuum. Section ^ concludes with some remarks concerning the renormalization process. 

There have also been a number of attempts to calculate some of the required renormalization factors 
numcr ically H 0. The renormalization factors computed from numerical simulations should agree with 
the perturbative calculations as the weak coupling limit is approached. This is an important test of the 
numerical renormalization techniques, which has not yet been attempted. The renormalization of the current 
has never been computed numerically. 



II Continuum HQET 

Heavy Quark Effective Theory is a way of studying a single heavy quark in a hadron when the mass of the 
quark is much larger than Aqcd- See Neubert [[IE) for a nice review of HQET. Mannel et al. make 
rigorous Georgi's |17| intuition that the heavy quarks at different velocities do not interact. They do so by 
showing that the QCD Green functions which involve two heavy quarks at different velocities go to zero in 
the infinite mass limit. So, there is a separate field for each heavy quark at each velocity. In HQET |iq] , the 
connection between the HQET fields and the quark fields in QCD (Q) is: 



h v (x) = e lMv - x P + Q{x) 
H v (x) = e lMv - x P-Q{x) 



Mvx-d (4) 



where P± = ± 1/). The new form of the QCD Lagrangian has h describing massless degrees of freedom 
and H describing fluctuations with twice the heavy quark mass. Further, explicit Gaussian integration of the 
H fields produces the effective, non-local Lagrangian. Upon integrating out the heavy degrees of freedom, 
the H term is replaced by a local term involving the light degrees of freedom, h, and the mass of the heavy 
quark, M. The Lagrangian is then expanded in the reciprocal of the heavy quark mass; the zeroth order 
HQET Lagrangian is 

£ c ff = h v iv ■ Dh v , (5) 

with the additional terms treated perturbatively as higher order in the reciprocal of the heavy quark mass. 
At zeroth order, i.e. in the infinite mass limit, the theory is independent of the mass of the heavy quark, and 
the Isgur-Wise function is universal (flavor-blind). 

In HQET, the momentum of the heavy quark (Mqv) is distinguished from the momentum of the light 
quarks and gluons (k, the "residual momentum"): 

Mv = Mqv + k. 

The residual momentum is the difference between the momentum of the hadron [Mv) and the momentum of 
the heavy quark. The velocity of the heavy quark becomes a parameter of the theory and it is the residual 
momentum which becomes conjugate to the position. In the infinite mass limit, the momentum of the hadron 
is due only to the heavy quark. 

The matrix element in the continuum M S scheme is connected to the matrix element calculated on the 
lattice by [§]| 

(v | V„ | v')^ = U<u | | vT t% = Zf(v | V, | V ') latt (6) 

where Z| is a continuum perturbative factor, is the lattice perturbative factor, and Z^ 1 is the ratio of the 
two. Falk et al. [fl9| calculated the continuum renormalization factor: 



Zl = 1 + {2 [1 - (v ■ v')r(v ■ »')] ln^/A) 2 + <S C } , (7) 



2 
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where 

\n(w + Vw 2 - 1) 
r{ - w > = / 2 = ( 8 ) 

and A is the gluon mass introduced as a infrared regulator. The dependence on A must cancel in Z£ , the 

ratio of Z^'s, of Eq. (|J). In the MS scheme, S c = [M. The calculation of the lattice renormalization 
factor, Z^, is the subject of the next section. Z^ will be discussed further in Sec. |^ when we discuss the 
matching from the lattice to the continuum. 



Ill Lattice HQET 

The Euclidean formulation of the lattice HQET action was introduced by Mandula and Ogilvie || . 

S = ^ Lotf(x) At^x) - i^v^ix) Aj + 2 A ~ j 4>(x)\ ■ (9) 

There is some freedom in the choice of which lattice derivatives are used in Eq. (^|). The tadpole improved 
finite differences are defined by 

& ip 3 = -3E±R _ ^ s (10) 

A_u^;j = x ' ipg-fi (11) 

u 



such that 



tpt t A ^ipgt is a forward difference, 

ipt t A_^ipg t is a backward difference, 

| (^A^t + ipl t A.^ St ^ is a centered difference, 

and Mo is the tadpole improvement factor [gjj. The tadpole renormalization of the lattice HQET action is 
subtle because of the constraint on the velocity; these subtleties are addressed in Appendix [X]. 

The centered difference approximates the continuum derivative to 0(a 2 ) (where a is the lattice spacing); 
both the forward and backward difference derivatives have O(a) corrections to the continuum. Therefore, it 
seems that the centered difference is the preferred type of derivative. This is true for the spatial derivative; 
however, Mandula and Ogilvie Jl3| emphasize that for consistency an asymmetric time difference must be 
employed, rather than a centered difference. If a centered difference is employed, then the propagator 
vanishes on alternate sites in the positive time direction and there is no continuum limit. The source of this 
problem is that the heavy quark fields are defined separately from the heavy antiquark fields and are distinct 
for each velocity [recall Eq. (Eh]; thus, heavy quarks can only propagate in one temporal direction. 

The lattice HQET action originally proposed by Mandula and Ogilvie [|| used a forward time derivative. 
The backward time derivative can be less convenient for use in simulations because a three dimensional matrix 
must be inverted for each time step. The forward time derivative only requires a matrix multiplication at 
each time step, and so is computationally cheaper to simulate. 

This choice of a forward time derivative has also been discussed by Davies and Thacker in the context 
of NRQCD. However, recent NRQCD calculations follow the prescription of Lepage et al. [^3| who use a 
backward time derivative but avoid having to invert a large spatial matrix by splitting the spatial part of the 
action over two adjacent time slices. Their action, which can be 0(a) improved, is symmetric with respect to 
time reversal, yet avoids the problems of the centered difference. Improved heavy- Wilson actions |E3 also go 
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over to the backward derivative in the static limit. Similarly, better choices for the HQET lattice actions can 
be made, and if we were to rerun the program with higher-order corrections, it would indeed be advantageous 
to use the backwards time derivative as is done for heavy- Wilson and modern NRQCD actions. But for our 
present purposes, the zeroth-order action suffices, and at this order the forward difference provides a technical 
advantage in computation. 

Since Aglietti's |l2j perturbative calculation used the backward-difference time derivative, we do the 
perturbative calculations for both types of time derivative. We can check our results against Aglietti's, 
against the results from the static theory p5[ , and also the static limit of NRQCD Comparison in 

perturbation theory between the forward and backward difference actions for the static case has led to 



the introduction of the "reduced wavefunction renormalization" discussed in Sec. [II. D and summarized in 
Appendix |B]. (Please see Appendix |C] for a comparison of the notation between the groups.) 

We introduce the notation 

-1 forward difference , 10 * 
-1 backward difference 

in order to compare the forward versus backward difference actions. Both the forward and backward differ- 
ence actions can be represented simultaneously by replacing A t by A CTt , where A CTt is either a forward time 
difference or a backward time difference, depending on the choice of action. 

Feynman rules can be derived from the action: 



quark 
propagator 

gluon 
propagator 

vertex 



tadpole 
vertex 

internal 
integrations 



v a 



A(fc) 



— e lapi - 1^ +V ^sin(p,) 



-I -1 



v 



<W (a 9 ^{T a ) bd (T a ) dc e^A ( $-^(T a ) bd (T a ) dc S m Pj 

\ I u ) t -r i \ 2 u Q 

*' a d 4 k 



(13) 

(14) 

(15) 

(16) 
(17) 



The T a are the color generators and Cf = 3 is the Casimir invariant. A is a gluon mass, which is needed to 
regulate the infrared divergences (as is done in the continuum) and which will be taken to zero at the end 
of the calculation. 

From the Feynman rules, it is straightforward to derive the usual self-energy tadpole [E tad (p)], 

and vertex [V(p,p')] corrections [the self-energy is £(p) + S tad (p)]: 



^tad 



(p) 



9 2 Cf 



2 2p j +k j 
2 



(2tt) 4 A(-fc) 




v cj (jL e to(P*+*d - l) + Y,j ^ sinfe + 



E^sinfe 



(2tt) 4 A(k) 



(18) 



(19) 
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V(0,0) = 



i 2 C f 



d 4 k 



(2tt) 4 A(-fc) 



V V c i<7k 4 /2 c iak 4 /2 

ul 



V C7 



- 1 



« 



\- v 3 v j k 3 kj 
> 7T- cos — cos — > / 



U 



(20) 



It is sufficient for our purposes to evaluate the vertex function with zero external momentum. The explicit 
p dependence is kept in the self-energy since the derivative will be considered. The integral which appears 
in the tadpole correction is standard and has the value p^yr / d 4: k/A(k) — 0.154933. 

The evaluation of the integrals is non-standard because of the problems caused by the spectrum of the 
Euclidean HQET action not being bounded from below. We follow the formalism developed by Aglietti 
and by Mandula and Ogilvie Jl3| ], in which we must first perform the k^ integration analytically and do 
so by transforming to z-space (z = e ±lki ). A contour is chosen that enforces the forward propagation of 
the HQET quarks [n3| as described below. (The connection to Minkowski space via a Wick contraction is 
discussed by Aglietti]l2j.) The resulting three dimensional integrals are then calculated numerically. (All 
of the numeric integrations were computed with the VEGAS routine |^| .) 

The analytic ki integration of Eqs. ( ^8[ ) through ( p0| ) reduces the four dimensional integration to a three 
dimensional integration. It is, however, more convenient to do this as a contour integration in z-space |l3| 
after an action-dependent change of variables 



For this change of variables, the gluon propagator is written 

A(fc M ) = 2^(1 - cos fcju) + {a\) 2 = 2 (1 - cos/c 4 ) + A 3 (fc) 



(21) 



(22) 



which defines As(k). 

The &4 contour (along the real axis) transforms into the unit circle in complex z-space. A subtlety arises 
when deciding which poles to enclose by the contour. The quark propagator pole appears as 



ZQ 



— p -wp4 



sinfe 



(23) 



The gluon propagator poles appear at 



A 3 1 
z± = 1 + — ± - 
± 2 2 



4A, 



(24) 



where A3 is defined by Eq. (22). The contour separates the gluon poles. The contour should enclose the 
quark pole and one of the gluon poles. The subtlety is in choosing which gluon pole. Because the energy- 
momentum relation from the quark propagator, Eq. (B4) of Appendix [b|, can be negative, we split fc-space 
(or z-space) into a positive-energy region and a negative-energy region and enclose the gluon pole which lies 
in the positive energy region of the space. For negligible external momentum with a quark momentum p + k, 
the upper k^ half-plane is positive energy and, using Eq. ( p"2| ) to distinguish the actions, it is convenient to 
define z via Eq. ( |2~lj ) such that z = e +lki for the forward difference action and z = e~ lki for the backward 
difference action. (For p— k, the lower k^ half-plane is positive energy and it is convenient to use z = e~~ l<Tki .) 
With either of these choices, the backward difference action will have the positive energy region outside of the 
z-space unit circle and the forward difference action will have the positive energy region inside the z-space 
unit circle. 
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The quark pole, 



Z Q 



1 - a- 



v ■ k 



(25) 



with positive-energy [using Eq. (B4), zq w 1 — ae] is just inside (outside) of the unit circle for the forward 
(backward) difference action. Since y/ A| + 4A3 > A3, we find z+ outside (and z_ inside) of the unit 
circle. z a (which is equal to z + for the forward difference and z_ for the backward difference) is therefore 
in the negative energy region. Since the k + (z+) gluon pole is always in the positive energy region for the 
backward-difference action and always in the negative energy region for the forward-difference action, we 
can write 



Backward Difference 
Forward Difference 



Z- 
Z- 



positive energy pole 
negative energy pole 
negative energy pole 
positive energy pole 



Z-r, 



negative energy pole 
positive energy pole 



(26) 



In both cases, it is the quark and the positive energy gluon poles which are enclosed by the contour regardless 
of where the quark pole actually appears. When the quark pole moves into the negative-energy region, it is 
necessary to deform the contour to keep the quark pole enclosed. (This is discussed for fc-space by Aglietti [Q 
and for z-space by Mandula and Ogilvie |l3j.) However, to simplify, one can equate this to the negative 
of the contour integral which encloses only the negative energy gluon pole. The three-dimensional integrals 
resulting from the contour integration have an action-dependent form due to the appearance of the negative 
energy gluon pole (z a ). This pole is a function of k. 

In order to compute the renormalizations, the unrenormalized propagator is compared to the renormalized 
propagator. (We include the mass term in order to calculate the mass renormalization.) The renormalized 
propagator has the form 

ifP(fc) = 7 T . (27) 



iu^k^ 



3 v)k, + AP + 0(k 2 ) 



The renormalization factors are obtained by Taylor series expanding the unrenormalized propagator: 



iH(k) = 



We used 



v a — e 

Uq 



1 1 + y — sm (%) + m o - E(fc, v) 



E(fc, v) = E(0, v) + fc 4 Xt + k j x j + °( fc2 ) 



(28) 



(29) 



and 



where 



Uq 



1 + iokn — Iiiuq + 0(k 2 ) 



dY(k) 



dk„ 



k=0 



(30) 



(31) 



Equation (|3C|) was used for the static case |27j, |28| to elucidate that the tadpole factor uq results in mass 
renormalization rather than wavefunction renormalization. Notice that uq has the perturbative expansion 
[1 — 7r 2 + 0[.g 4 )], so In wo ~ 0{g 2 ); the higher order terms are neglected. After a little algebra which 
involves the addition and subtraction of some deducible terms, one can write the propagator in the form 



iH(k) = {[I- 5Z] 



i (vq + Sv Q ) fc 4 + 



E 

3 



it 



kj + (M + 6M) + 0(k 2 ) 



(32) 
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which implies the expressions for the renormalizations: 



SM = 
SZ = 



Sv 



Z- 1 

r U * 

U 

Uq 

v r o - v 



— E(0, v) — crwo l n wo 
-roo^4 — Uo w j Aj 

-i«o— A 4 - (1 + vf)Xi - Vi VjXj 



Uq 



-i(vo - l)Xi - v Q u ^ VjXj 



(33) 
(34) 

(35) 

(36) 



We make the following points regarding these expressions: First, in the HQET formalism, the residual 
momentum is conjugate to the position, leaving the velocity as a free parameter. As discussed by Aglietti , 
the velocity is renormalized on the lattice. In the continuum, the four- vector X^ is proportional to v^, the 
only available four- vector; this implies that there is no velocity renormalization. On the lattice, with reduced 
rotational symmetry, this is not the case. Secondly, if uq is set to unity and the special case of v — v z z is 
taken, then these reduce, for the backward-difference case, to Aglietti's result JO). Thirdly, 6^5 is a notation 
to remind the reader that this quantity renormalizes — rather than Vj as can be seen in Eq. (p5[). For Uq = 1, 



our St^ corresponds to Aglietti's Svj. Further, the velocity renormalization can be written as follows: 



r,tad 



6, tad rytad 



r,tad 6, tad r^tad 

V Q — V Q £ VQ 



ztad 




z 



t ad 



8vq 
v 



(37) 



g 2 c F , 

16ir 2 



It is taken 



Finally, the uq that appears in these expressions is the perturbatively expanded uq = 1 

at lowest order (unity) and the terms higher order in g 2 are ignored because X^ ~ C(g 2 )- The result is that 
the wavefunction renormalization and the first term of the mass renormalization, S(0, u), are the same to 
0(g 2 ) whether or not one uses tadpole improvement. The velocity renormalization is affected by tadpole 
improvement as 



Z 



t ad 



s— 



9 2 C f 



1 + ir-T6^^) 



UO 

,2 



9 2 C F 
16tt 2 



[c(v) 



(38) 



where the (— n 2 ) is from the perturbative expansion of uq and c(v) is the same [to C{g 2 )] regardless of tadpole 
improvement. 

Of the renormalization factors [Eqs. ( |33| ) to (|3^)], only the mass renormalization, Eq. (33), depends explic- 
itly on the choice of forward or backward time difference (the a parameter). However, all the renormalization 
factors implicitly depend on a via the X^ functions. The explicit dependence of the mass renormalization on 
a is zero when tadpole improvement is not used; this is discussed further in Appendix [B], above Eq. (B1C). 



III. A Velocity Renormalization 

Mandula and Ogilvie p"3] ] renormalize Vj = 
will be renormalizing v, and calculating c(v) defined by Eqs. (|3 



rather than v. We will not be using their notation, rather we 



and IXii 



16tt 2 



c(vj). 



(39) 
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This parallels the notation of Aglietti [ fL2| . (Sec Appendix |c] for a comparison.) Recall that this is the 
perturbative renormalization to the tadpole-improved velocity. Neither Mandula and Ogilvie nor Aglietti 
use a tadpole-improved action. 

The expression for c(v) is found from the self-energy Feynman diagrams as expressed through Eq. (|35|). 
Continuing to use the a — ±1 to distinguish between the actions, we find 



civ) = — 



d 3 k (-^z a {k) + i + vf) sm(k t ) + £ sin(%) 



7T 7 V /A i + 4A 3 



a(^(fc)-i) + E,| sin (^) 



<?( k ) ~ J2j ^ cos 2 (^) z CT (fc) - + «? ) cos(fc») - £ i?4i «f cos(fc,-) 



(40) 



The uq are perturbatively expanded such that at this order in g 2 , they can be replaced with unity. (They 
are included as a reminder that in the next order there will be an effect.) Note that z a {k) is the negative 
energy gluon pole, defined by Eqs. (^4j) and (|2^), introduced from the residue of the contour. 

Mandula and Ogilvie p3[ perform an expansion in small velocity and present the velocity renormaliza- 
tion as coefficients to powers of the velocity. (This is convenient in that whenever a calculation at a new 
velocity is desired, the value for the velocity has precalculated coefficients so that the calculation need not 
be done repeatedly.) While this is straightforward for the velocity renormalization, the divergences in the 
wavefunction renormalization and the vertex correction make this technique more complicated for these other 
calculations. However, if we consider the expansion for the velocity renormalization, then we get consistent 
results at C(v 6 ) (notice that our format is slightly different because c renormalizes v rather than v) 



c(v) 



= C 00 + C 2 oavf + CQ20 

+0400^ + C 2 2(A 2 ^ + c 040 



Y «* + c «22 Y J2 vfrl 



~6 ~4 \ ~2 ~2 \ T ~4 ~2 \ T \ ^ ~2 ~2 

-CGOOVi +C 4 20«i Z^ V 3 + C 240Wi 2^, V 3 + C 222«i }^ Z^ V j V k 



+C060 + c 042 Y Y + ■ ■ 



(41) 



The forward and backward difference results of this expansion are listed in Table [l]. Mandula's and Ogilvie's 
results are reproduced by the first two columns. Our results for the same special case (backward difference, 
v x = v y = 0) that Aglietti considers [jl2| are listed in Table and agree with Aglietti where they overlap. The 
three columns of the forward difference are: c(v) according to Eq. (|4C|), its expansion through sixth-order 
in small velocity according to Eq. (fll]), and its expansion through second-order of the velocity expansion 
using only the first three terms of Eq. (|4f|). The latter confirms Mandula's and Ogilvie's result; however, 
the sixth order result (using the coefficients of Table [T]) is in much better agreement with the exact result 
(as one would expect). Although Table ||only considers motion along a single axis, our more general results 
indicate that for the forward difference action it is sufficient to use the velocity expansion to sixth order. 

For the more general case of all the spatial velocities not equal to zero, we present the results for the 
forward difference action at small velocities in Table ^. This is the factor, c z {v), which renormalizes the z 



component of the velocity according to Eq. 



The renormalizations for the v x and v y components can be 



deduced from the table by symmetry. Notice that the v z renormalization is affected by each component of v, 
not merely by v z . The numerical size of the perturbative factors in Tables ^ and ^| are both large. Tadpole 



improving the perturbative factors, by adding 
size of the perturbative contribution. 



-7r to them as in Eq. (Bq) , does not substantially reduce the 
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Table 1: The coefficients, c mn /, used in the velocity renormalization when expanded in powers of the 
velocity to C(v e ) according to Eq. (fi"l]). s is the order of the velocity term, found by summing the indices: 
s = m + n + l. The first set is for the forward difference action; the second set is for the backward difference 
action. If you consider the velocity in only one direction, then only the top row is relevant. 





Forward Difference 




Backward Difference 




n I 


s =0 s=2 s=4 


s—6 


s=0 s=2 s=4 


s=6 





-28.07(3) -4.977(6) -1.093(3) 


-0.458(2) 


11.78(1) 0.33(2) -0.88(3) 


-2.03(3) 


2 


-4.292(6) -2.100(6) 


-1.380(6) 


10.26(2) 9.49(6) 


7.0(2) 


4 


-1.010(3) 


-1.346(6) 


7.62(3) 


28.1(2) 


2 2 


-1.005(6) 


-1.36(1) 


9.53(6) 


43.4(3) 


6 




-0.469(2) 




39.98(6) 


4 2 




-4.54(2) 




109.6(6) 



Table 2: This table lists the velocity renormalization for both forward difference (our choice) and backward 
difference actions for the special case v x = v y = 0. The last two columns solve the expanded equation 
through the superscripted order. The c{v) entries are exact, that is, not expanded in the velocity. Note that 
the v — > 0.0 limit is considered even though there is no need to calculate the renormalization coefficient when 
v = 0. (c has no interpretation in the static limit.) 







(Backward) 




(Forward) 








c(v) 


c{v) 


c (6) 


c (2) 


c(v 


-> 0.0) 


11.779(4) 


-28.06(1) 


-28.06(1) 




c(v 


= 0.1) 


11.899(5) 


-28.40(1) 


-28.40(1) 


-28.38(1) 


c(v 


= 0.2) 


12.275(5) 


-29.44(1) 


-29.44(1) 


-29.42(1) 


c(v 


= 0.3) 


12.966(5) 


-31.35(1) 


-31.35(1) 


-31.29(1) 


c(v 


= 0.4) 


14.036(7) 


-34.39(1) 


-34.39(1) 


-34.28(1) 


c(v 


= 0.5) 


15.67(1) 


-39.17(1) 


-39.17(1) 


-38.95(2) 


c(v 


= 0.6) 


18.05(1) 


-46.90(1) 


-46.90(1) 


-46.44(2) 


c(v 


= 0.7) 


20.82(3) 


-60.44(2) 


-60.45(2) 


-59.47(3) 
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Table 3: The velocity renormalization, c z (v), for the forward difference action for several general (small) 
velocities. The uncertainty is at most 2 in the last digit. 





0.00 


0.05 


0.10 


0.25 


0.00 


-28.06 


-28.15 


-28.40 


-30.14 


0.05 


-28.16 


-28.23 


-28.47 


-30.19 


0.10 


-28.39 


-28.47 


-28.72 


-30.46 


0.25 


-30.12 


-30.20 


-30.45 


-32.17 




0.00 


0.05 


0.10 


0.25 


0.00 


-28.16 


-28.23 


-28.48 


-30.19 


0.05 


-28.26 


-28.34 


-28.55 


-30.27 


0.10 


-28.49 


-28.56 


-28.82 


-30.53 


0.25 


-30.20 


-30.30 


-30.54 


-32.26 


V x \ 


0.00 


0.05 


0.10 


0.25 


0.00 


-28.40 


-28.48 


-28.72 


-30.45 


0.05 


-28.49 


-28.55 


-28.81 


-30.53 


0.10 


-28.72 


-28.80 


-29.06 


-30.78 


0.25 


-30.45 


-30.52 


-30.79 


-32.50 


e.V» 


0.00 


0.05 


0.10 


0.25 


0.00 


-30.15 


-30.21 


-30.46 


-32.20 


0.05 


-30.22 


-30.31 


-30.56 


-32.27 


0.10 


-30.48 


-30.54 


-30.80 


-32.54 


0.25 


-32.21 


-32.27 


-32.51 


-34.26 
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To give an idea about the magnitude of the velocity renormalization, we consider (3 = 6.0 with \v\ = 0.5, 
and use the bare lattice coupling. The non-tadpole improved multiplicative factor Z v , = 0.67; the cor- 
responding tadpole improved number is Z^ = 0.75. If the boosted coupling, <? 2 / u o ^l), is used then 
>gtad _ q 59 s i p e f the Isgur-Wise function essentially depends quadratically on the velocity renor- 

malization, this makes perturbation theory unreliable to analyze the simulation data and thus numerical 
renormalization techniques must be used j|, ||, [l3|, |l4| . 



III.A.l Aside: Slow HQET 

In Aglietti's (l^] initial calculations, the velocity renormalization was presented as a function of the velocity. 
However, Mandula and Ogilvie fl3|| expanded the velocity renormalization in a power series in the velocity, 
which allowed them to compare their perturbative results with the numbers from their numerical renormal- 
ization technique. The expansion of the renormalization factor in velocity seems to be similar to Aglietti's p<| 
idea of slow HQET, where the v ■ D term is a perturbation on the static theory. Slow HQET was studied in 
perturbation theory by Aglietti and Gimenez p0| , where they demonstrated that slow HQET agreed with 
HQET in the infrared and ultraviolet limits. It would be interesting to understand the connection between 
slow HQET and the HQET formalism of Mandula and Ogilvie. 

We have found expressions for the velocity renormalization in terms of the coefficients for the backward 
difference action (Table [l]) and note that the C042 coefficient of the backward difference is rather large, at 
109.6(6), - much larger than the equivalent coefficient for the forward difference action. This could indicate 
a problem with the expansion for the backward difference; the forward difference coefficients, which we 
checked through 0(v 6 ), are all reasonably close to unity. Aglietti and Gimenez do not calculate the sixth 
order coefficient for the velocity renormalization (although they take the other renormalizations to this order); 
however, the renormalization factors are only required to quadratic order in the velocity in order to compute 
the slope of the Isgur-Wise function from simulations of slow HQET. (The slow HQET formalism is used to 
directly calculate the derivatives of the Isgur-Wise function, using the "moments" technique [^9| . Aglietti 
and Gimenez found that the expressions for the higher order derivatives of the Isgur-Wise function, 
beyond the slope, contained operators that diverged with an inverse power of the lattice spacing and that 
must be subtracted off in the simulation.) 



III.B Mass Renormalization 

As with the velocity renormalization, we define x(v) as the mass renormalization without the <7 2 -prefactor. 
For comparison, Aglietti [[l2| also does this; however, we prefer (for comparison to the static limit of the 
forward difference NRQCD theory) to have our x{v) proportional to +S(0, v). So, our x is the negative 
of Aglietti's. (See Appendix ^ for a comparison between groups.) We also include the effect of tadpole 
improvement. 

SM = — 2(0, v) — (tvq ln uo = — — — — av o lnuo (42) 

I67H a 

Recall that lnu is C(g 2 ). From Eq. (||), 

m = 2 »°/ , £t . -^> + s4 c ° 82 (») . ..^ ,„„,„■,,, (43) 

« J y/E$ + 4A 3 [„ (J.,„ ( t) - 1) + Ej Jt.,i„(fc,) 

where the 871" 2 (0.154933) is from the tadpole contribution (E tad ) which is partially canceled by the second 
term in Eq. ( |42| ) as it should be. The Uq are again perturbatively expanded such that, at this order, they 
can be replaced with unity. They are included here as a reminder that in the next order they will have an 
effect. The values of this integral are listed in Table ^ in the next subsection. As they are relevant to the 
reduced wavefunction, wc will discuss these there. 
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III.C Wavefunction Renormalization 

The results of the wavefunction and reduced wavefunction renormalization can be compared not only to 
Aglietti @ (backward difference, HQET) and Mandula and Ogilvie £jj| (forward difference, HQET), but 
also to Eichten and Hill 25, [uj (backward difference, static theory) and the static limit of Davies and 
Thacker (forward difference, NRQCD). (Appendix [c] compares the notations between groups.) 

Recall that the wavefunction renormalization can be found as Eq. (|34|). During this calculation, as with 
the velocity renormalization, the ko integration is done analytically with the same comments as were made 
earlier. This introduces the residue from the negative energy gluon pole, z a {k). Again using the a = ±1 to 
distinguish between the actions, the result for SZ is 



5Z 



g 2 c F H 



d 3 k 



-2a^(fc)+E,|sin(%) 



16tH tt J ^Ai + 4A 3 [ [a^ik) - 1) + £ sm (fc 



a(±z a (k)-l) + Z/t si ^ k 



16tt 2 



8tt 2 (0.154933) 



where again the 87r 2 (0. 154933) is from the tadpole contribution and the uq are again perturbatively expanded 
such that at this order, they can be replaced with unity. This three-dimensional integral has a logarithmic 
divergence. The way with which this is typically dealt is to add and subtract an integral with the same 
logarithmic divergence which is solvable analytically. We call this integral SZ C and use the small k limit 
because we are interested in the infrared (low energy) divergence. The difference SZ — SZ C is finite and 
calculated numerically. SZ C (found analytically) will have a finite piece, which is added back to the numerical 
calculation, as well as a divergent piece. The divergent piece contributes to the coefficient of the ln(A 2 a 2 ) 
term in the renormalization of the lattice Isgur-Wise function. 

Although the "continuum-like" limit of SZ is actually 



g 2 C F 2v 2 f d 3 k 
16tt 2 ~ J 2Vfc 2 + A 2 c 



"0 ^ 



+ 



(VP + A 2 a 2 ) + E, Jfc, 



+ 



9 2 C F 
16tt 2 



8tt 2 (0.154933), 



the first and third terms are finite. Since we are interested in the infrared divergent piece, we will define SZ C 
as the second term. By taking advantage of uo| P ert = 1 + 0(g 2 ) as well as by using the velocity normalization 



we get 
SZ C - 



g 2 C F 1 [ R d 3 k 
~^^vlJ Vk 2 + X 2 c 



2 r 



The upper limit R is arbitrary because this term is added and subtracted. Interestingly, this is the same 



Ill LATTICE HQET 



14 



integral for both actions. The result of this integral is (with v = v'j. 



9 2 C F 1 
16tt 2 vl 



A<Cfi 



In 



9 2 C F 
16tt 2 



In I 



Vi? 2 + A 2 a 2 -R v{l~v 2 ) \ Vi? 2 + A 2 



Vi? 2 + A 2 a 2 + iR 



vR 









[ 21n U0 


-?ln( 











(44) 



The divergent piece is —2 ln(A 2 a 2 ). This is the wavefunction contribution to the divergence in the renormal- 
ization of the Isgur-Wise function. We can, for convenience, set R = 1%. 

As with the mass renormalization, the results for the wavefunction renormalization are listed in Table ^, 
but discussed further in the next subsection. Note that, as with x(v) for the mass renormalization, the 
wavefunction renormalization is referred to by e(v) and defined by 



SZ 



9 2 C F 
16tt 2 



-21n(AV) +e(v)] 



(45) 



32 



on 



III.D Renormalization of the Reduced Wavefunction 

The perturbative factors for various heavy quark effective field theories depend subtly |l^. ^2], 
whether the forward or backward time derivative is used in the action. It is expected that in the Euclidean 
formulation, the propagator as a function of time and the residual 3-momentum (i.e. Fourier transforming 
&4 into t) will have the dependence e~ £t — e~ mt . However, it turns out (App. |b|) to have the dependence 
A e -m(t-a) w }j ere fj — ±i distinguishes the actions. Eichten and Hill [0 noticed this relation and found 
that if one fits, instead, to A'e~ mt (where A' = Ae m<y ) this changes the wavefunction renormalization by 
subtracting (or adding) the mass renormalization. It also "reduces" the wavefunction renormalization to a 
common answer for both the forward and backward difference actions. Since it is convenient to fit to e~ mt 
and the reduced value is the same for both actions, this is a popular choice. Unfortunately, in lattice HQET 
away from the static limit the reduced values (for the forward and the backward-difference cases) are not 
equal, as we will show. 



Equation ( p2[ ) defines x(v) in terms of the mass renormalization. Equation 
the wavefunction renormalization. Appendix H derives Eq. (Bll), 



e (v) = e(v) — ax(v)/vQ e 



3) defines e{v) in terms of 



(46) 



which is the relation for the reduced value of the wavefunction renormalization. The tadpole term is in 
x(v) [it gets canceled in the mass renormalization of Eq. (|42)] and, as noticed for the static case in |2?], p|, 
the wavefunction and reduced wavefunction renormalizations remain unaffected by tadpole improvement. 
Table f| lists our values for these functions for a large velocity range. Notice that in the static limit, the 
reduced value for the two actions is the same. This is the expected result. Also notice that our forward 
difference value for e(0.0) agrees with Davies and Thacker |2^] (in their notation C = Z + aA = —15.4). Our 
backward difference table agrees with Aglietti fll^l ; and the static limit of the backward difference action, 
e(0.0) = 24.44, is also in agreement with Eichten and Hill pq| . While it is still convenient to use the reduced 
result and fit to e~ mt in the static limit, the forward and backward difference will have different reduced 
wavefunction renormalizations away from the static limit. 



III.E Vertex Correction 



The vertex correction also has differences between the actions and a divergence which must be subtracted as 
was done for the wavefunction renormalization. However, this has the further complication that it depends 
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on the velocities of both the incoming and the outgoing quarks. So whereas the wavefunction renormalization 
is a function of v, the vertex correction, SV(v,v'), is a function of the initial and final velocities. 

After analytically doing the contour integration over the k± variable and dealing with the poles as dis- 
cussed previously, we find 



5V(v,v') 



g 2 C F 2 f d 3 k 
16tt 2 n J y/Al + 4A 3 



v&vo , 
u 



v cr \—z a (k) -1 ] + V — sin(fcj) 
\u J u a 



u' a (-z a {k) sm(kj) \ -(47) 



The Uq are once again perturbatively expanded such that at this order, they can be replaced with unity. 
(They are included here as a reminder that in the next order they will have an effect.) For the rest of this 
subsection, we will explicitly set uo=l. 

As claimed by Aglietti [Q , this must have the form 



g 2 C F 2 

1671" 2 7T 



[2 (vv')r(vv') ln(AV) +d(v,v')] 



(48) 



The lattice coefficient must have this form if it is to cancel the continuum divergence lnC/x) (which was 
computed by Falk et al. |fl9|). Of primary interest is that it be a function of v ■ v', the only non-trivial 
invariant constructible from the heavy quark velocities v and v' . We find that numerically the lattice 
divergence coefficient agrees with the continuum divergence coefficient — these are listed in Table |. 

Aglietti only gives results for the v = v z z with v' — case. We have found a problem with an equation 
which he uses and have introduced a better expression for finding the divergence. We have also extended 
the calculation to the forward difference action and to nonzero v' . 



III.E.l Vertex Correction with v' 







In the v' — > limit, we can reproduce Aglietti's results; the numerical values are listed in the next subsection. 
Originally, we could not satisfactorily reproduce Aglietti's numbers, especially for the v = 0.0 result. As 
we investigated this, we found a problem with the A — > limit, specifically there was a subtlety with the 
interchange of limits (A — > versus v — * 0). We believe that Aglietti's choice of integral subtraction can be 
improved. This subsection discusses this subtlety. Tables [| fj], and || were produced with our choice. 

Since the SV integral is divergent, a technique similar to that used for SZ can be used; however, it needs 
to be modified because the "continuum- like" limit of this integral, 6V C , is not analytically manageable. 
However, a second integral, 6V CC , can be taken such that SV — 5V C and 5V C — 6V CC are each finite. These 
numeric integrals are then done separately and added together along with the finite piece of 5V CC . 

Aglietti refers to 5V C as /, and our SV CC is analogous to his L. Aglietti uses the notation 5(v) = L — I; 
we will use the analogous definition S'(v) — 5V CC — SV C . To be explicit, in the small-fc limit Eq. ( |4T|) reduces 
to 



6V c (v,v') 
5V c (v,0) 



d 3 k 



g 2 C F 2 

16tt 2 tt J 2Vfc 2 + A 2 <r J 

g 2 C F 2 f d 3 k 



Vfc 2 + A 2 a 2 + V - Vjkj Vfc 2 + A 2 a 2 + V • v'fy 



16tt 2 tt J 2(fc 2 + A 2 a 2 ) 



Vfc 2 + A 2 a 2 + vjk. 



(49) 



(50) 
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Table 4: Mass [#(£>)] and wavefunction [e(v)] renormalization functions for the backward-difference (bd) 
and forward-difference (fd) actions. The bd numbers reproduce Aglietti's table The bd and fd numbers 
for x(v) and e'(v) should agree only in the static limit (v=0). e'(v) = e(v) — ax(v)/vo are the reduced 
wavefunction. Notice that we define x(v) as the negative of that of Aglietti. In addition, v x = v y = 0. 





Backward Difference 


Forward Difference 


V 


x(v) 


e(v) 


e'(v) 


x(v) 


e(v) 


e'(v) 


0.0 


-19.92(3) 


24.43(4) 


4.53(1) 


-19.93(1) 


-15.40(1) 


4.530(4) 


0.1 


-19.87(3) 


24.64(4) 


4.875(4) 


-19.99(1) 


-15.75(1) 


4.141(2) 


0.2 


-19.69(3) 


25.24(4) 


5.97(1) 


-20.17(1) 


-16.82(1) 


2.935(2) 


0.3 


-19.34(3) 


26.36(4) 


7.91(1) 


-20.47(1) 


-18.78(1) 


0.759(4) 


0.4 


-18.75(3) 


28.14(4) 


10.96(1) 


-20.97(1) 


-21.91(1) 


-2.694(6) 


0.5 


-17.72(3) 


30.94(5) 


15.60(2) 


-21.72(2) 


-26.83(2) 


-8.015(8) 


0.6 


-15.79(3) 


35.44(5) 


22.82(4) 


-22.89(1) 


-34.74(2) 


-16.44(2) 


0.7 


-11.15(3) 


44.2(1) 


36.27(10) 


-24.79(2) 


-48.56(4) 


-30.85(2) 



Table 5: The coefficient of the lattice divergent ln(Aa) piece must and does reproduce the continuum 
divergent coefficient [4 (v ■ v') r(v ■ v')] in order to correctly cancel the In A. Errors are at most three in the 
last digit shown. In addition, v x = v y = v' x — v' = 0. 



-t\ v 

v' \ 


0.00 


0.10 


0.20 


0.30 


0.40 


0.50 


0.60 


0.70 


0.00 


4.001 


4.014 


4.050 


4.127 


4.237 


4.396 


4.618 


4.953 


0.10 


4.012 


4.004 


4.015 


4.059 


4.139 


4.269 


4.462 


4.753 


0.20 


4.053 


4.011 


4.002 


4.018 


4.064 


4.160 


4.319 


4.573 


0.30 


4.126 


4.058 


4.014 


4.000 


4.016 


4.079 


4.194 


4.412 


0.40 


4.238 


4.137 


4.066 


4.016 


4.000 


4.021 


4.095 


4.257 


0.50 


4.398 


4.267 


4.159 


4.077 


4.022 


4.000 


4.026 


4.133 


0.60 


4.622 


4.460 


4.316 


4.197 


4.095 


4.028 


3.998 


4.040 


0.70 


4.956 


4.753 


4.571 


4.412 


4.257 


4.135 


4.041 


4.000 
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Aglietti approximates SV c (v,v' = 0) with 
g 2 C F 2 f d 3 k 



L 



-1 



16tt 2 tt J 2(k 2 + X 2 a 2 ) 



g 2 C F -2^1 

167T 2 V 



1 



A 2 a 2 



(51) 



However, this gives a A-dependent 5{v). In spherical coordinates, S(v) has the form 

9 2 C F 2 



167T 2 v J k 2 + X 2 a 2 



hi I 



VA: 2 + A 2 a 2 + vk 
y/k 2 + A 2 a 2 — ufe 



hi 



(52) 



So long as v is finite, we can take the limit as A — ► 0. However, if we want both v — > and A — > 0, a problem 
arises: the result in the limit A — * is not the result at A = 0. This is a case in which the limits cannot be 
interchanged. To be rigorous, we break up the integration into a region for k < A and for k > A. 



(5(6 — > 0, A — > 0) 



16?r 2 ' 



lim 



dk 



k' 2 



dk 



\{k 2 + X 2 a 2 ) 
k 2 k 



{k 2 + X 2 a 2 ) 



4 lim 



(k 2 + \ 2 a 2 f 2 {k 2 + X 2 a 2 ) 
1 



i+ H 4 a& 



g 2 C F 

16tt 2 a^o 
+ limO( A2a2 / i? 2)+ lim 0{%) 

A^0 A>e^0 



4 h HA^ 



(53) 



(54) 
(55) 



While a A divergence was expected for SV, the difference S(v) must be finite. 8(v — > 0, A — » 0) is infinite 
because the logarithms do not cancel exactly. It happens that (5(A) has a minimum around A w 10 -5 ; 
at this value, if v is taken to zero, then Aglietti's d(v = Q,v' = 0) = —4.53 can be calculated from 
lims^o [6V(v,0)-6V c (v,0)+ finite part of V cc (i,0)} = -5.75 and S(i -> 0, A « 10" 5 ) = -1.22. How- 
ever, Eq. ( |54| ) clearly shows that (5(A — > 0) blows up. This can be seen for A < 10~ 5 . 

To avoid this problem, we write Eq. ( |50|) as 



JV c ({i,w' = 0) 



g 2 C F 2 

167T 2 7T 



d 3 fc 



-1 



1 



f ; ■ k 



2(k 2 + X 2 a 2 f 2 

which expands as follows and allows a better definition of 5V cc (v, v' = 0): 

g 2 C F 2 /" d 3 k -I 



(56) 



5y c (w,w' = o) 

SV cc (v,v' = 0) 



16tt 2 tt J 2(k 2 + X 2 a 2 f 2 
g 2 C F 2 f d 3 k 



167T 2 7T 



-1 



2(fc 2 + A 2 a 2 ) 



1 + 5 • k 



(57) 
(58) 



We find that this makes 6'(v) stable to small A and that it generalizes to give useful results when v' ^ 0. In 
the equation analogous to Eq. (p4j), the logarithms cancel and the result is finite in the A — ► limit. 



III.E.2 Vertex Correction with v' ^ 

This case requires the continuum- like expression for Eq. ( [47| ) . Recall Eq. (|4^) . Again there are problems if we 
use Aglietti's trick of setting A to zero in the factors with v and v'. The problems are: (1) the A-dependence 
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Table 6: The finite piece of the backward difference vertex correction d(v, v') for v x = v y = v' x = v y 



v' \ 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.0 


-4.528(2) 


-4.579(2) 


-4.756(2) 


-5.088(3) 


-5.642(3) 


-6.592(4) 


-8.454(5) 


0.1 


-4.583(2) 


-4.459(2) 


-4.446(2) 


-4.567(3) 


-4.892(3) 


-5.562(3) 


-7.018(4) 


0.2 


-4.755(2) 


-4.447(2) 


-4.228(3) 


-4.131(3) 


-4.186(3) 


-4.523(3) 


-5.525(4) 


0.3 


-5.088(3) 


-4.570(3) 


-4.126(3) 


-3.768(3) 


-3.514(3) 


-3.457(4) 


-3.890(4) 


0.4 


-5.640(3) 


-4.890(3) 


-4.183(3) 


-3.517(3) 


-2.882(4) 


-2.320(4) 


-2.007(5) 


0.5 


-6.598(4) 


-5.556(3) 


-4.523(4) 


-3.460(4) 


-2.320(4) 


-1.071(5) 


0.283(6) 


0.6 


-8.452(4) 


-7.022(4) 


-5.529(5) 


-3.884(4) 


-2.003(5) 


0.280(6) 


3.322(8) 


)le 7: 


The finite piece of the forward difference vertex correction d(v, 


v ) lor v x = v y 


/ / 

= V X = Uy 


v' \ 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.0 


-4.528(2) 


-4.513(2) 


-4.471(2) 


-4.402(3) 


-4.280(3) 


-4.102(3) 


-3.822(4) 


0.1 


-4.516(2) 


-4.527(2) 


-4.509(2) 


-4.462(3) 


-4.370(3) 


-4.221(3) 


-3.971(4) 


0.2 


-4.471(2) 


-4.508(2) 


-4.513(3) 


-4.487(3) 


-4.437(3) 


-4.315(3) 


-4.108(4) 


0.3 


-4.402(3) 


-4.457(3) 


-4.494(3) 


-4.505(3) 


-4.472(3) 


-4.396(3) 


-4.224(4) 


0.4 


-4.283(3) 


-4.371(3) 


-4.431(3) 


-4.474(3) 


-4.481(3) 


-4.446(4) 


-4.325(4) 


0.5 


-4.108(4) 


-4.219(3) 


-4.316(3) 


-4.398(4) 


-4.443(4) 


-4.459(4) 


-4.402(4) 


0.6 


-3.822(4) 


-3.973(4) 


-4.112(4) 


-4.226(4) 


-4.320(4) 


-4.398(4) 


-4.424(4) 



is incorrect (which implies that the difference is A-dependent), (2) the limit as v' — > does not reproduce the 
results of the previous subsection, and (3) the integral is too difficult. So, once again, we will try to retain 
the A-dependence as follows: we approximate 



5V c (v,v') 



9 2 C F 2 

167T 2 7T 



d 3 k 



by 



5V CC («,«') 



2(fc 2 + A 2 a 2 ) 3/2 



d 3 k 



v ■ k 



Vk 2 +\ 2 a 2 



v' -k 



Vk 2 +\ 2 a' 2 



-i + E 



g 2 C F 2 

J 2 (fc 2 + A 2 a 2 ) 3/2 h + « . k] h + V . k 



While this does solve both the A-dependence problem and the ? -t problem, it only barely solves the 
difficulty of the integral. However, in spherical coordinates, it allows the |fc| integral to be separated from 
the angular integration. We can solve this integral by doing the \k\ integration analytically. (Since this is 
where the A-divergence exists, it is the only piece that needs to be done analytically anyway.) Having thus 
removed the ln(A) term, we can numerically calculate and add back the angular integration along with the 
finite piece of the |fc| integration. 

Since this is symmetric in v and v' , the results for the vertex correction should be also. Our extension to 
nonzero v' shows that this is the case: Tables ^ and ^ show our results for backward and forward difference 
vertex correction at general velocities. Notice that the results are symmetric about the diagonal, v — v' . 
Notice also that the first row and the first column of Table |6| both reproduce the backward-difference results 
of v' = in Table ||. Table [?] shows the results for the forward-difference vertex correction at general velocities 
and the first row and column reproduce the forward-difference results of Table [|. 
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IV Lattice to Continuum Matching 

To renormalize the Isgur-Wise function, which is proportional to the current in Eq. (|^), we need the current 
renormalization, which can be assembled from the wavefunction and vertex renormalizations calculated in 
the previous section. This involves 



Z%(v)Z v (v,v')Z%(v') 



l + \5Z Q (v) 



1 + SV(v, v') 



l + \sz Q (v') 



(59) 



So, following Aglietti's lead |T^], we define 

f(v,v') 
f'(v,v') 



l + -[SZ Q (v)+8Z Q (v')]+SV(v,v') 



\[e{v)+e{v')]+d{v,v') 
\[e'{v)+e'{v')]+d{v,v') 



where a reduced Isgur-Wise correction, /', is defined using the reduced wavefunction, e', which was used 
with a fit model of the form e~ mt . Since the wavefunction reduction does not affect the vertex correction, 
d, the perturbative factor for the lattice Isgur-Wise function is: 



Z \ = l + T7T^ {~ 2 V) r(v ■ v')} ln(Aa) 2 + f(v, «')} 



12tt 



(60) 



where the divergences have been isolated to calculate the finite pieces and r(v ■ v') has been defined by 
Eq. (||). If we did not wish to use the reduced value, the divergence would stay the same and we would 
merely replace the /' with /. We have already shown (Table ||) that the divergent piece of the lattice vertex 
correction cancels exactly with that of the continuum; thus the lattice logarithm coefficient is written with 
the same form as for the continuum correction, Eq. (]?]). Now the continuum correction, Z|, can be divided 
by the lattice correction, Z 1 ^, to find the lattice to continuum matching factor 



Zf(v, v f ) = l + JL-{2[l-(v> v>) r(v ■ v')} ln(^a) 2 - f(v, v')} 



(61) 



The expression in Eq. (|6l]) is suitable for renormalizing the Isgur-Wise function extracted by taking ratios of 
two and three point functions Q. However, to improve statistics, HQET simulations extract the Isgur-Wise 
function using ratios of three point functions only ||, [li], |i| . We discuss this additional complication below. 

Our results for d, f, and /' are listed in the following tables. Recall Tables and ^ show our results 
for backward and forward difference vertex correction at general velocities. Table lists our results for the 
vertex correction, the current correction, and the reduced current correction in the backward and forward 
difference actions for v' = 0. The backward difference reproduces Aglietti's results. Tables |] and |l0| show 
our results for backward difference current and reduced current corrections at general velocities. Again, the 
results are symmetric about the diagonal. Tables [H] and ^ for the forward difference current and reduced 
current corrections at general velocities are also symmetric about the diagonal. Notice that the first rows 
and columns of Tables ^, [?], and ^ through [l^ reproduce Table ||. Notice also that although the different 
actions give the same result in the static limit (v — > 0, v' — > 0), this is not the case at any other velocity. 

For continuum HQET in the MS renormalization scheme at zero recoil (v ■ v' = 1), Z'i = 1 and the finite 
piece is zero [[ll|. This corresponds to the diagonal (v = v') of the tables which contain our results. On the 
lattice, however, if the conserved current is not used, f'{v,v) is not constrained to be zero. We account for 
this next. 

To deal with the finite piece of the renormalization, we note that the numeric extraction of the Isgur- 
Wise function on the lattice does not calculate the Isgur-Wise function directly. The numerical extraction is 
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Table 8: The finite piece of the vertex correction, d(v,v'), of the forward- and backward-difference actions 
for v x = v y = and v' — 0. The backward-difference action results should and do reproduce Aglietti's table 
up to the correction mentioned in the text of this paper. Also listed is the current correction f(v,v') = 
\e{v) + 2 e {v') + d(v, v') and the reduced current correction f'(v, v') = \e'{v) + ^e'(v r ) + d(v, v') which form 
the correction for the lattice Isgur-Wise function. 



Backward Difference Forward Difference 



v d(v,v' = 0) f(v,v' = 0) f'(v,v' = 0) d(v,v' = 0) f(v,v' = 0) f'(v,v' = 0) 



0.0 


-4.526(2) 


19.92(1) 


0.000(2) 


-4.527(2) 


-19.94(1) 


0.000(2) 


0.1 


-4.578(2) 


19.96(1) 


0.122(2) 


-4.511(2) 


-20.09(1) 


-0.174(2) 


0.2 


-4.757(2) 


20.08(1) 


0.489(2) 


-4.474(2) 


-20.58(1) 


-0.740(2) 


0.3 


-5.089(2) 


20.33(2) 


1.129(2) 


-4.401(2) 


-21.50(1) 


-1.755(2) 


0.4 


-5.639(4) 


20.63(2) 


2.100(4) 


-4.282(2) 


-22.93(1) 


-3.364(2) 


0.5 


-6.597(4) 


21.09(2) 


3.459(4) 


-4.104(4) 


-25.21(1) 


-5.844(2) 


0.6 


-8.432(6) 


21.50(2) 


5.23(1) 


-3.800(4) 


-28.89(1) 


-9.755(4) 


0.7 


-14.80(1) 


19.53(4) 


5.57(2) 


-3.354(6) 


-35.32(2) 


-16.529(8) 



Table 9: The finite piece of the backward difference current correction, f(v,v'), for v x = v y = and 



v' \ 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.0 


19.903(8) 


19.958(8) 


20.100(9) 


20.324(10) 


20.654(9) 


21.093(9) 


21.51(1) 


0.1 


19.961(8) 


20.195(10) 


20.485(8) 


20.921(10) 


21.516(10) 


22.245(10) 


23.03(1) 


0.2 


20.087(8) 


20.498(9) 


21.022(9) 


21.662(9) 


22.522(9) 


23.565(10) 


24.84(1) 


0.3 


20.320(9) 


20.947(9) 


21.681(9) 


22.57(1) 


23.73(1) 


25.20(1) 


27.02(1) 


0.4 


20.646(10) 


21.494(9) 


22.495(9) 


23.74(1) 


25.25(1) 


27.24(1) 


29.78(1) 


0.5 


21.090(9) 


22.228(9) 


23.579(10) 


25.175(10) 


27.20(1) 


29.86(1) 


33.48(1) 


0.6 


21.499(10) 


23.017(10) 


24.84(1) 


27.02(1) 


29.75(1) 


33.47(1) 


38.75(1) 



Table 10: The finite piece of the backward difference reduced current correction, f'(v, v'), for v x = v y = 
and v' = v' = 0. 



A' 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.0 


0.004(3) 


0.118(2) 


0.491(2) 


1.134(3) 


2.098(3) 


3.458(4) 


5.228(6) 


0.1 


0.119(2) 


0.422(3) 


0.975(2) 


1.821(3) 


3.021(3) 


4.675(4) 


6.822(6) 


0.2 


0.488(2) 


0.974(2) 


1.733(3) 


2.810(3) 


4.273(4) 


6.258(5) 


8.858(6) 


0.3 


1.128(3) 


1.822(3) 


2.812(3) 


4.142(4) 


5.922(4) 


8.292(5) 


11.484(7) 


0.4 


2.096(3) 


3.026(4) 


4.274(4) 


5.922(4) 


8.073(5) 


10.954(6) 


14.875(7) 


0.5 


3.459(4) 


4.672(4) 


6.250(4) 


8.291(5) 


10.957(6) 


14.504(7) 


19.477(9) 


0.6 


5.219(6) 


6.839(6) 


8.856(6) 


11.477(7) 


14.893(7) 


19.470(9) 


26.10(1) 
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Table 11: The negative of the finite piece of the forward difference current correction, —f(v,v'), for 
v x = v v = and v' = v' = 0. 



5'\ 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.0 


19.923(6) 


20.086(6) 


20.591(7) 


21.500(7) 


22.953(7) 


25.223(8) 


28.897(9) 


0.1 


20.095(7) 


20.282(7) 


20.786(7) 


21.714(7) 


23.206(7) 


25.513(8) 


29.217(8) 


0.2 


20.591(7) 


20.784(7) 


21.342(7) 


22.297(7) 


23.810(7) 


26.136(8) 


29.882(10) 


0.3 


21.491(7) 


21.742(7) 


22.298(7) 


23.268(8) 


24.819(8) 


27.188(9) 


30.987(10) 


0.4 


22.941(8) 


23.192(8) 


23.794(7) 


24.824(8) 


26.394(9) 


28.829(9) 


32.66(1) 


0.5 


25.220(8) 


25.498(9) 


26.138(8) 


27.177(8) 


28.804(9) 


31.283(10) 


35.15(1) 


0.6 


28.885(9) 


29.22(1) 


29.897(9) 


30.993(10) 


32.64(1) 


35.16(1) 


39.16(1) 


0.7 


35.32(1) 


35.68(1) 


36.41(1) 


37.59(1) 


39.30(1) 


41.90(1) 


46.01(2) 



Tabic 12: The negative of the finite piece of the forward difference reduced current correction, —f'(v,v'), 
for v x = v v = and v' = v' = 0. 



v' \ 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.0 


-0.004(3) 


0.178(2) 


0.737(2) 


1.755(2) 


3.365(3) 


5.843(3) 


9.773(4) 


0.1 


0.181(2) 


0.380(2) 


0.968(2) 


2.008(2) 


3.645(3) 


6.154(3) 


10.113(4) 


0.2 


0.740(2) 


0.967(2) 


1.583(2) 


2.649(3) 


4.311(3) 


6.849(3) 


10.851(4) 


0.3 


1.759(2) 


2.009(3) 


2.652(3) 


3.742(3) 


5.441(3) 


8.020(4) 


12.061(5) 


0.4 


3.366(3) 


3.642(3) 


4.312(3) 


5.442(3) 


7.173(3) 


9.796(4) 


13.881(5) 


0.5 


5.841(3) 


6.154(3) 


6.858(3) 


8.022(4) 


9.796(4) 


12.450(5) 


16.618(6) 


0.6 


9.772(4) 


10.118(4) 


10.855(4) 


12.059(5) 


13.886(5) 


16.610(6) 


20.845(8) 


0.7 


16.525(6) 


16.911(6) 


17.687(6) 


18.964(7) 


20.839(8) 


23.675(9) 


27.996(10) 
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more manageable using the technique of Mandula and Ogilvie J34[ where the ratio of the three-point quark 
propagator, G (defined explicitly in our concurrent numerical paper ||), gives a ratio of lattice Isgur-Wise 
functions: 



Av v' G v ' v '(t)G v '' v (t) 



(v q +v^G^(t)G-'^(t) \Z(v,v)\ 
This technique exploits the continuum normalization of the Isgur-Wise function at zero recoil 

Z(vv)=Z(l) = l. (63) 



(62) 



Since v 2 is normalized to 1, the denominator of Eq. ( |62[ ) can be set to unity in the continuum. This ratio 
also allows the normalizations and smearing-function dependence to cancel, so we expect that 



Zf{v,v')£} att («,«') 



Zf(v,v)£ latt (v,v) 



Zf(v',v)£ latt (v',v) 



Zf(v',v')C att (v',v') 



\Z cont (v -v')\' 

|£ COnt (l)| 2 



= \e° nt (v -V')\ 



Thus, our unrenormalized calculation of 

£ratio («,«') 

must be renormalized by 



>y) = 



£ latt (u,y)£ latt (y',y') 



' Zf{v,v')Zf{v',v) 
Zf(v,v)Zf(v',v') 



written as 



Zf atio (v, «') UUv, v') r nt (v ■ v>). 



(64) 
(65) 



On the lattice, £} att {v,v) does not obey Eq. (g3|) unless a conserved current is used; nevertheless, 
6atio («,«') (by definition) acts like the continuum Isgur-Wise function even if the conserved current is 
not used. Without the conserved current, £} att (v,v) ^ 1, but the normalization cancels in the ratio so that 
£ratio(w, v) = 1. Thus, Z^ atio will be symmetric in v and v' and will have the property Z^l tio (v, v) = 1. 

Expanding Eq. (p4), we find 



7 cl — i 

^ratio — 1 



lg 2 C F 



{2 [(1 - (v ■ v') r(v ■ v')) + (1 - {v 1 ■ v) r(v' ■ v)) 



2 167T 2 

— (1 — (v • v) r(y ■ v j) — (1 — (v 1 ■ v') r(y' ■ v'))] ln(yua) 2 
-f'(v,v') -f'(v',v) +f(v,v) + /'(«',«')}• 



Using v ■ v — v' ■ v' = r(l) = 1 and f'(v, v') = f'(v', v), this reduces to 



Zf atio (v,v') = l + 



12tt 2 



2 (1 — (v ■ v') r(v ■ v')) ln(/ia) 2 — f'(v, v) + 



f'(v,v) + f(v',v') 



(66) 



(67) 



which not only has the correct divergent coefficient but we also see a new finite piece which is manifestly 
zero on the diagonal. The wavefunction renormalization cancels explicitly in Eq. (67), so /' can be replaced 
by the vertex correction d. 
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V Conclusions 

We have calculated the renormalization of the lattice b — > c current by considering the lattice Isgur-Wisc 
function. This calculation extends previous work by including tadpole improvement, by extending to non- 
zero initial and final velocities, and by considering forward as well as backward difference actions. 

By considering the forward difference action and the backward difference action side-by-side, we find 
non-trivial differences between the two. The practical difference in a lattice calculation is that the backward 
difference requires a matrix inversion at each step of the calculation. The differences in the renormalization 
are that the gluon poles over which one integrates are interchanged; away from the static limit, the reduced 
values are no longer equal; and the velocity renormalization, when expanded as powers of the velocity, stay 
small for the forward difference, but grow large for the backward difference. 

Of greater concern is that the velocity renormalization is not terribly small. We have shown that the 
velocity renormalization can be expanded in small velocity and that the coefficients remain on the order 
of unity at higher orders (at least for the forward difference action). These coefficients are given here to 
C{v 6 ). The nonperturbative calculations are giving smaller renormalizations |l3|, [l4| and these should be, in 
principle, more reliable. This should be considered in more detail, especially the slow HQET for the forward 
difference action. 

Although our results confirm other groups' calculations where they overlap, the integrals and divergences 
are subtle and must be managed with care. When we combine our renormalizations into a current correction 
with the ratio introduced by Mandula and Ogilvie |34|] , such that the finite piece of the current correction is 
—f'(v, v') + 3 [f'(v, v) + /'(«', £')]> we nn d that all of our results have the appropriate limits and cancelations. 
These expressions are used in our concurrent numerical paper [g| to compute the slope of the Isgur-Wise 
function using lattice HQET. 



A Appendix: Tadpole Improvement 

Tadpole improvement is a mean field improvement which (at lowest order) cancels the effects of the 
large "tadpole" Feynman diagrams. In the HQET, there is no coefficient (analogous to k in the Wilson 
action) which is common to both Ut and Uj and which allows one to a posteriori tadpole improve any 
previous calculation which was not tadpole improved. Fortunately, as noticed by Mandula and Ogilvie |l3[| , 
the evolution equation can be written such that the uq is grouped with ij — Vj / Vo . Thus, tadpole-improved 
(tad) Monte-Carlo data can be constructed from the non-tadpole-improved (nt) data by replacing v nt — > u tad 
and by including two overall multiplicative factors ( v o/ v ± & d was not included by Mandula and Ogilvie): 

7 ,nt 

G tad M tad ,^ ad ) = Uo-'^M"^) (Al) 

Vq 

In addition, the tadpole-improvement of a simulation requires adjusting the velocity (analogous to adjusting 
k) according to w tad = uov nt , subject to the normalization ( w tad ) 2 = i and (v nt ) 2 = 1. The adjustment on 
the velocity is then 

«r = tf[i+(i-ui)(vf)']- lh (A2) 

vf d = u vf[l + {l-ul){vf) 2 ]-^. 1 ' 

The tadpole improved data is at a velocity which is shifted from the original tadpole unimproved data. 
Previous HQET calculations have either not included tadpole improvement [jl2| or have had difficulties with 
it |l3| . Although one should start with a tadpole improved action, we find it convenient to be able to tadpole 
improve a calculation a posteriori because there are choices for how one can determine the mean-field value 
uo PI. 
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B Appendix: Reduced Renormalizations 



One can define a "reduced" wavefunction renormalization and relate it to the fit-model exponential. We 
begin by considering the propagator as a function of time t and the residual momentum fcj^] 



iH(t,k) 



dk,A 



a ik 4 



2" 



9 



Since iff ~ e , the energy- momentum relation can be found: 



e = —aln 



1 — <7 — sin(fco) 



— sinffc,) 



(B3) 



(B4) 



Aglietti |12[ notes that the energy goes to zero for both k = and k = w, but provides a physical argument 
for why this doubling problem has a negligible effect in the HQET. 



In Eq. (B3), it may be noticed that the 0-function has a different argument for the different actions. 
Though it was phrased differently, this was also noticed by Davies and Thacker p3] who give recursive 
expressions for the Green function evolution equation for the two cases of a forward or a backward difference 
in their NRQCD action. 

In order to consider the renormalization effects of the fitting form, consider the next loop-order of the 
propagator as a function of the time and the residual 3-momentum, 



dk 



2tt 



|ii?(fc 4) fc) +iH(k 4 ,k) S(fc) iH(k 4 ,k)\ 



(B5) 



Following Aglietti [[[2) , we will make use of 

E(fc) = £(0) + £4X4 + ^ fe,-Xj + 0(fc 2 ) 

j 

= -5M tad + [-av a (l-5Z) + 5v ]\n(u ) 



-SZ 



v (T — e 

Uq 



-SvqCT I — e 
u 



E 

j 



Vj 

6— } sm(kj) 



(B6) 



where SM is the tadpole improved mass renormalization (versus SM nt the not tadpole improved mass 
renormalization) defined by 



5M tad = -£(0,5) - crv lnu = SM nt - crw kiu . 



(B7) 



It may also be noticed that since lnuo ~ Oig 2 ), the [(vqcSZ + 8vo)\nuo] can be neglected as G{g A ). We 
further note that terms of the residual momentum, C{k), can be neglected.n (The residual momentum can 
be adjusted by introducing a "residual mass." ) Finally, we note that the 5v and 6^ can be collected with 

1 Recall that the residual momentum, rather than the full momentum, is conjugate to the position. 
2 The calculation including these terms is available from author JC. 
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the bare velocity in precisely the proportion necessary to renormalize each velocity. To solve these integrals, 



one needs to put Eq. (B5) into a form which allows the use of 

dx e %ax 



-co 27r ( e<X 

dx e l 



, 2ir (e lx - l) 2 
With these relationships, we find (eventually^) 



= 9(a) (B8) 
= (o-l)e(o). (B9) 



zH^(t,k) = e(t + 1 ^-) U °f + SZ) exp{-(t - a) [M g ]} (l + Oy 9 2 ) + Oy 2 )) 



(l + ^) 



where M a is an action-dependent function of the renormalizations, of the velocity, and of the momentum; 
and vq (l + Sv °/v ) = ^o™- The relevant point is that, as was said previously, for the forward difference action 
one should fit to a form of exp(— Mf[t — 1]); whereas for the backward difference action one should fit to a 
form of exp(— Mj,[t + 1]) [i.e. fit to exp{— M G [t — a)}]. However, if one chooses to fit to the form exp(— Mt), 
then the coefficient Z = (1 + SZ) gets changed to Ze aM ° w Z(l + aM a ) « (1 + 5Z + oM a ). To C{g 2 ), 
neglecting Oik) terms, M„ = (SM tad + av r cn lnu )/v r cn = -£(0)/^ cn [Recall Eq. @]. So, to this order, 
the "reduced" wavefunction renormalization is 



(BIO) 



Z' = Z-a^{Q)/vl cn 

= (1 + 5Z - aS(0)/^ en ) 

= [ ~ 2 ln(AV) + e{i) - ax ^y< n ^ ) ■ 

This is also written in terms of the finite pieces 

e'(v) = e(v) - ax{i)/v™ n . (Bll) 

The tadpole term is in x(v) [it gets canceled in the mass renormalization of Eq. (|42])1 and, as noticed for 
the static case in [ p7| pq| , the wavefunction and reduced wavefunction renormalizations remain unaffected 
by tadpole improvement. 



C Appendix: Notation 

When comparing between the results of HQET, NRQCD and the static theory, the difference in notation 
starts to become a factor. Where Davies and Thacker (NRQCD) used A for S(0), Aglictti (HQET) uses 
A(p) for the non-tadpole portion of the self-energy as well as using A for a particular grouping of terms for 
convenience in the calculation. We are going to maintain Davies' and Thacker's use of A and give new names 
to Aglietti's As. However, since Aglietti considers the velocity-dependence of various quantities, we will use 
Aglietti's notation for a variety of velocity-dependent functions. The velocity will be relevant for the HQET, 
but not for the static theory nor for the NRQCD. In the HQET, the functional dependence is on v defined 
by 

- = E^ 2 = E5- (012) 

J 3 

Note that Aglietti calls this u. 

3 We found that there are two © terms. One goes as 0^t + ~ — j = 0i, the other as Q^t + ~ — ^ = 83. We resolved this 
assuming we were interested in late enough times (t > —1) that 63 = 9\ = 1. 
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Table 13: Comparison of notation between Aglietti JL2J and Mandula and Ogilvie p3| . Note also that 
Aglietti only considers motion in the z-direction. Finally note that in the last row, Mandula and Ogilvie 
consider 8vi, but Aglietti considers Sv z (u, — Vi — Vi / V0 )- To convert between the two, one must include a 
factor of v 2 . 

Mandula and Ogilvie Aglietti Comparison 

X — -iX 4 X X\ Ag = X 4 \ MO = iX \ MO 

X i Y Y \a 9 = X 3\mO 

Vi v z . _ . 

vi = — Uz = ~ u Aa 9 = v Amo 

v o v o 

Svi = -—(Xi-ViXo) Sv z = -ivQV z X-v\Y &v_ _ 

= -^(Xj+ivjXi) = -v 2 (Y + iu z X) ~ ~ V °~ 



Aglietti calls the mass renormalization SM; he also puts in a negative sign, which we leave out. Aglietti 
notes that for the HQET, this is velocity dependent, and defines a function a(u) which is proportional to his 
SM 

g 2 C F s(«)Iab 



5M\j 



Vq does not appear in NRQCD and is 1 in the static limit. 



In calculating the wavefunction renormalization, 



dp,! 



(C13) 



is needed. Mandula and Ogilvie use the 



p=0 



notation X^ . This is a useful notation and does not conflict with either Davies and Thacker or Eichten and 



X \ MO and Y\ Ag 



X: 



See Table [13] for an explicit comparison. 



Hill. Aglietti names these as A^ 
We choose to use Mandula's and Ogilvie's notation 

In the definition of the velocity renormalization, there is a further subtlety. Mandula and Ogilvie consider 
Svi = x^™' — Vi, but Aglietti considers Sv z = ui ron ' ) — v z (with the definition Ui = Uj = Vi / Vo )- As shown in 
Table [l3|, a factor of Vq must be included to translate between ^ anc i §£. . J n addition, Mandula and Ogilvie 

include the pre- factor \Fz in their definition of c(v) in Eq. (|39|). 

Mandula and Ogilvie do not calculate the wavefunction renormalization, therefore we will compare Agli- 
etti's wavefunction renormalization to Davies and Thacker (while using Mandula's and Ogilvie's notation for 
Xfj). Aglietti uses SZ = Z — 1 for the wavefunction renormalization. To relate this to Davies and Thacker, 
we note that 



Z \a 9 = 1 



Xo-5>*i = 1 



9 2 C F 
16tt 2 



!ln(oA) +e(v) 



(C14) 



where the ln(aA) term comes from doing the self-energy integral. It is 5Z = Z — 1 which is Davies' and 
Thacker's C. 



16tt 2 



21n(aA) 2 + e(v) = C\ DT = Z\ DT + aA\ 



I DT 



(C15) 



In addition, because of some discrepancies discussed in Sec. pII.D , it will be convenient to define a "reduced 
value of e," [e^(v) = e'(v)]: 



C F 
16tt 2 



-21n(aA) +e'(v) = Z\ 



IDT 



(C16) 



This reduced value can be found from e(v) and x(v) as expressed in Eq. (Bll) 
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